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Abstract 

We show that the problem of solving recurrence relations for L-loop 
{R + l)-point Feynman integrals within the method of integration by 
parts is equivalent to the corresponding problem for [L + i?)-loop vac- 
uum or {L + R— l)-loop propagator-type integrals. Using this property 
we solve recurrence relations for two-loop massless vertex diagrams, 
with arbitrary numerators and integer powers of propagators in the 
case when two legs are on the light cone, by reducing the problem to 
the well-known solution of the corresponding recurrence relations for 
massless three-loop propagator diagrams with specific boundary condi- 
tions. 
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1 Introduction 



The method of integration by parts (IBP) [|T], 0] within dimensional regularization |^ 
is based on integrating by parts in dimensionally regularized Feynman integrals in 
loop momenta and putting to zero surface terms. Although this procedure has been 
justified only for off-shell Feynman integrals |Q, no examples are known where it 
breaks down so that it is successfully applied to Feynman integrals on the mass shell, 
at threshold and to various integrals in the loop momenta that arise in asymptotic 
expansions of 'standard' Feynman integrals in limits of momenta and masses. 

The idea of the method is to use equations following from IBP in order to de- 
rive recurrence relations that reduce any Feynman integral from a given family to 
a set of master integrals, e.g. with lowest powers of propagators. If this program 
is fulfilled then the evaluation of the given family of integrals reduces to algebraical 
manipulations. 

In this paper, we show that different problems of solving IBP relations are in 
fact equivalent so that one can use known solutions to provide solutions for other 
families of Feynman integrals. We show that the problem of solving IBP relations for 
L-loop {R + l)-point Feynman integrals is equivalent to the corresponding problem 
for (L + _R)-loop vacuum or {L + R — l)-loop propagator-type integrals. 

We apply this property to reduce solution of the IBP relations for two-loop mass- 
less vertex diagrams, with arbitrary numerators and integer powers of propagators 
in the case when two legs are on the light cone, to the well-known solution for 
massless three-loop propagator-type diagrams, but with specific boundary conditions. 

2 Equivalence of recurrence procedures 

Our goal is to prove the equivalence (in some special sense) of the IBP relations 
for multi-loop Feynman integrals with the same number of independent kinematical 
invariants. To do this let us present (following p, with updated notation) the IBP 
relations in a modified form. We start with L-loop vacuum Feynman integrals with 
the maximum possible number of the propagators, = L{L + l)/2: 

where pi {i = 1, . . . , L) are loop momenta and Da = A^Jpi ■ pj — ^a^n^ (a = 1, . . . , A^; 
summation over repeated indices is understood). We use dimensional regularization 
with (i = 4 — 2e. The matrix A^^ is supposed to be symmetrical with respect 
to upper indices. This maximal number A^ of the propagators provides the possi- 
bility to express any scalar product of the loop momenta as a linear combination of 
the factors in the denominator. Feynman integrals for specific graphs are obtained 
for appropriate choices of the N x N matrices A]^ . Diagrams of practical interest 



1 



which usually have less number of the propagators can be considered as special cases 
with non-positive powers of some propagators. If the number of the propagators is 
greater than N, partial fractioning can be used to deal with integrals with at most 
propagators. 

IBP relations are obtained H| by acting by the operator (d/dpi) ■ pk on the 
integrand, putting the integral of the derivative to zero and expressing all the terms 
resulted from the differentiation through the initial family of the integrals. Starting 
from (|I]) we obtain 

d6ikf{ni,...,nN;d) = 2i^'(I; + . . . , n^v; c?) , (2) 

where .., n^, ...) = /(. . .,na-l, . . .) and . . .^.) = . . ,^^ + 1, . . .) 

(no summation in in the last relation ). The elements A^^ of the matrix inverse 
with respect to Af, i.e. with v4^' = + 6ii6jk)/2, arise when expressing the 

scalar products pi -pj in the numerator, which appear as a result of the differentiation, 
in terms of the denominators Da'- 

Pk-pi = A'J{D, + fi,). (3) 

Suppose now that the integrals depend in addition on R external momenta pi 
{L < i < L + R) and we are interested in their values at a particular point Pi ■ pk = 
fj.ikm'^- Without loss of generality we may assume that det(yUjfc) 7^ 0. Indeed, if this 
determinant is equal to zero then the momenta Pi at this point are linear dependent. 
Let the rank of the matrix fin. be R' < R. This means that we can represent all 
expressions as functions of only R' external momenta so that we come back to the 
(slightly simplified) original problem. 

The number of independent propagators (i.e. factors in the denominator) which 
can be constructed from L loop and R (independent) external momenta is A^^i = 
L{L + l)/2 + LR, and the number of additional ('external') invariants is N2 = R{R + 
l)/2. Let us expand the integrals in formal series (i.e. expand the integrands in 
the corresponding Taylor series) in the external kinematical invariants Pk ■ Pi, k,l = 
L + 1, . . . , L + R, OT, in other words, in the 'denominator-like' objects Da, a = Ni + 
1, .., Ni + N2, depending only on the external momenta: 

F{ni,...,nN;d;m'^;pL+i,...,PL+R) = J 

E (m2)^'^/2-^"-+^V(ni, 

In particular, the 'on-shell' value at Pi ■ Pk = fJ^ik^"^ of the integral is represented by 
the first term /(ni, .., tzat^, 1, .., 1; d) in the expansion in Dik = Pi ■ Pk — fJ^ikfn'^, while 
the other terms in the expansion play an auxiliary role when proving the equivalence. 

Acting by (d/dpi) ■ pk, {i = 1, . . . ,L;k = 1, . . . , L + R) on the integrand in 
(0) we obtain A^^i IBP relations exactly in the 'vacuum' form (0). To control the 



/■ dV...dV 

'"J D1\..DI'' 

N1+N2 

..,n^,+^,;d) n ^a""'- (4) 

a=Ni+l 
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evolution of the A^2 'external' indices we need additional N2 relations. We obtain 
them by acting by pk ■ (d/dpi), {i,k = L + 1, . . . , L + R) on (Q). These new rela- 
tions look like (Q) with the only exception that they have no term on the left-hand 
side proportional to the space-time dimension d. To transform these additional rela- 
tions into a manifestly equivalent form we rescale the initial integrals by a power of 
the determinant of external kinematics invariants normalized on its 'on-shell' value 
D = det{pi ■ pk)/ det(/ijfcm^) , L < i, k < L + R: 

F{m, . . . , n^; rf; . . .) = . . . ,n^; d; . . .) , (5) 

with F expanded in / by the same equation (H) as in the case of F. For pure 'on- 
shell' values, the prefactor in the right-hand side of (|^) is equal to one. For higher 
terms of the expansion in the external invariants, this prefactor is expanded and 
provides a linear substitution (easily invertible) of coefficients / throught /. As a 
result, the coefficients /(ni, . . . ,nAr^+7V2; d) exactly satisfy the 'vacuum' relations (||) 
with N = Ni + N2. 

Of course, boundary conditions for the recurrence procedures can be different. In 
particular, when solving the recurrence relations for pure the 'on-shell' L-loop 
point integrals we can use the combinatorial results for the corresponding (L+i?)-loop 
vacuum or {L + R — l)-loop propagator recurrence relations, but with the additional 
condition that integrals with non-positive values of the 'external' indices should be 
put to zero. 

3 Recurrence procedure for massless 2-loop vertex 
diagrams 

Let us apply the equivalence property described in the previous section to the eval- 
uation of the massless two-loop vertex diagrams shown in Fig. |I|a,b,c with two legs 




(a) (b) (c) 

Figure 1: Three types of two-loop vertex diagrams with general numerators and 
integer powers of propagators: (a) planar, (b) non-planar and (c) non-Abelian. 
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on the light cone, Pi = P2 = 0, general polynomials in the numerator and arbi- 
trary integer indices of the propagators. The diagrams depend homogeneously on 
= (^1+^2)^ = —Q^- For example, a general Feynman integral for the planar graph 
Fig. I^a can be written as 



{k ■ 0"^ 



(P - 2pi ■ A;)"3(A;2 + 2p2 ■ - /)2)"6 



(6) 



The standard iO prescription is implied, i.e. k"^ = k'^ + iO. We have chosen the 
irreducible numerator to be the scalar product of the two loop momenta, k ■ I. 

The planar and non-planar integrals with the powers of propagators equal to one 
and some numerators of low degree were evaluated in |^ in expansion in e. The 
planar integral with the powers of propagators equal to one was evaluated by IBP 
in gamma functions for general e in 0. Using IBP, an algorithm for the evaluation 
of general planar integrals without numerators has been developed when expanding 
two- loop vertex diagrams in the Sudakov limit [§. This solution of the recurrence 
relations was reproduced in 0. For planar diagrams with arbitrary numerators, a 
complete solution was presented in The non-Abelian diagrams Fig. |^c are non 
more complicated than the planar ones. The IBP relations for them are simple, and 
any diagram can be expressed in gamma functions for general values of d. Thus the 
only non-trivial case is in fact the non-planar diagram of Fig. |l]b. 

Following the results of Sect. 2 we reduce the evaluation of the diagrams of Fig. |l] 
to the evaluation of three-loop massless propagator diagrams of Fig. ^,b,c, with 
arbitrary numerators and integer powers of propagators. As is well known the latter 






(a) 



(b) 



(c) 



Figure 2: Three types of three-loop propagator diagrams with general numerators 
and integer powers of propagators: (a) planar, (b) non-planar and (c) Mercedes. 

can be reduced, using the algorithm of , to a linear combination of master integrals 
Fig. 1^ and recursively one-loop integrals Fig. ^ 

Using the statement justified in the previous section we conclude that to solve 
IBP relations for Fig. we can use recurrence relations for Fig. |^ presented in |^ 
in an explicitly solved form. The massless external legs should be identified with a 
pair of lines of the propagator diagrams which are connected to an (arbitrarily fixed) 
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Figure 3: Master diagrams for 3-loop propagators. 




Figure 4: Recursively one-loop diagrams for 3-loop propagators. Integer indices of 
lines and numerators are arbitrary. 

external vertex. The additional condition formulated above means that, when using 
the recurrence procedure, one should put to zero integrals with non-positive values 
of these 'former external' indices of the lines. 

To make a correspondence between the two given families of the diagrams let us 
choose, for definiteness, the right vertex in Fig. Ha,b,c and lines that are incident to 
this vertex so that the vertex diagrams of Fig. |l| will be obtained from the propagator 
diagrams of Fig. |^ by cutting the two right internal lines and, vice versa, the propa- 
gator diagrams will be obtained from the vertex diagrams by adding a new external 
vertex and a pair of lines which connect them with the two right external vertices in 

Fig. g. 

Thus, to evaluate a vertex diagram from Fig. ^ with a given set of indices n,, we 
can use the following prescriptions: 

• using the algorithm presented in |l| for the corresponding 3-loop propagator 
analog express it as a linear combination of the master integrals of Fig. |] and 
recursively one-loop diagrams Fig. ^; 

• omit the contribution of the recursively one-loop diagrams in the first raw of 
Fig.g 
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• substitute the values of the recursively one-loop diagrams from the second raw 
in Fig. I by the corresponding values of the recursively one-loop vertex-type 
integrals shown in Fig. ^ if the indices of the right two lines are equal to one. 
If at least one of the right indices is greater than one, then evaluate the prop- 
agator diagram and express the result (written in gamma functions) through 
the corresponding propagator diagram with both indices equal to one. Finally, 
perform the above substitution in this result; 

• substitute the values of the master integrals in Fig. |^ by the corresponding 
values of the master integrals shown in Fig. ^. 



Figure 5: Master diagrams for 2-loop vertices. All the lines have indices equal to 
one and the numerator is equal to one. 





Figure 6: Recursively one-loop diagrams for 2-loop vertices. Integer indices of lines 
and numerators are arbitrary. 



The correspondence between the two sets of the integrals is obvious: the diagrams 
from Fig. |^ and ^ are obtained respectively from Fig. |^ and ^ by cutting a pair of the 
right lines. Note however that the second master vertex integral in Fig. ^ is in fact 
recursively one-loop and easily evaluated in gamma functions for general e = (4 — (i)/2 
so that the only 'true' master vertex integral is the non-planar diagram in Fig. |^ which 
was evaluated in in expansion in e up to the finite part: 

d7 



{{k + /)2 - 2p^{k + IMk + ly + 2p2{k + - 2p^k){P + 2p2l)kn^ 

2 



1 TT^ 83C(3) 597rn {iT^'^'^e 



3e 



120 



{Q 



2\2+2e 
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where 7_b is the Euler constant. 
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